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Experimentalists  seeking  to  improve  the  coherent  lifetimes  of  quantum  bits  have  generally  focused  to  date 
on  addressing  their  respective  leading  decoherence  mechanisms  through,  for  example,  improvements  to  qubit 
designs,  materials,  and  system  isolation  from  extrinsic  perturbations.  In  the  case  of  the  phase  degree  of  freedom 
in  a  quantum  superposition,  however,  the  coherence  that  must  be  preserved  ultimately  includes  that  of  the  qubit 
relative  to  the  system  clock.  In  this  manuscript  we  clarify  the  impact  of  clock  instability  on  qubit  dephasing 
and  provide  quantitative  estimates  of  fidelity  upper-bounds  set  by  noisy  phase  fluctuations  in  the  clock.  We 
first  demonstrate  analytically  that  phase  fluctuations  in  the  clock  -  a  master  oscillator  typically  implemented  and 
referred  to  as  the  “local  oscillator”  (LO)  -  is  indistinguishable  from  a  pure  dephasing  field  arising  from  other 
environmental  mechanisms.  Using  these  results,  we  then  apply  commonly  quoted  LO  phase-noise  specifications 
to  rigorously  calculate  performance  bounds  on  qubit  operational  fidelities  due  to  LO  phase  fluctuations.  We 
find  that  for  state-of-the-art  performance  specifications,  phase  fluctuations  at  frequencies  far  from  the  carrier 
dominate  operational  error  rates.  Furthermore,  the  control  bandwidth  impacts  the  degree  to  which  thermal  noise 
forms  an  effective  phase-error  floor.  Our  results  motivate  additional  research  directions  needed  in  support  of 
quantum  information  systems  and  highlight  challenges  and  advantages  that  particular  qubit  technologies  possess 
in  the  context  of  phase-noise  induced  errors. 


I.  INTRODUCTION 

Dephasing  in  qubit  systems  is  commonly  induced  by  en¬ 
vironmental  fluctuations  of  a  qubit  bias  or  control  parameter, 
e.g.,  an  external  magnetic  field,  which  modulates  the  qubit- 
state  energy  splitting  and  hence  its  dynamic  evolution.  In  turn, 
experimentalists  have  accordingly  focused  their  efforts  on  re¬ 
ducing  both  the  level  of  environmental  noise  and  the  qubit 
sensitivity  to  that  which  remains.  As  a  result  of  these  efforts, 
the  performance  of  various  qubit  technologies  [1]  has  grown 
over  recent  years,  and  dephasing  times  are  now  measured  on 
timescales  approaching  milliseconds  or  even  seconds  [2-8] 
with  operational  fidelities  now  reaching  99.999%  [9-11].  At 
these  performance  levels,  noise  in  the  control  system  emerges 
as  an  important  contributor  to  the  error  rate,  and  it  will  eventu¬ 
ally  become  a  leading  dephasing  mechanism  as  the  quantum 
information  community  continues  to  mitigate  other  residual 
error  sources. 

In  this  manuscript  we  identify  dephasing  noise  induced  by 
master-clock  phase  fluctuations  as  a  potentially  significant 
source  of  error  in  today’s  qubit  systems  and  one  that  will  cer¬ 
tainly  become  more  prominent  as  error  rates  improve.  We 
express  the  effective  dephasing  Hamiltonian  in  terms  of  LO 
phase  fluctuations,  establishing  the  correspondence  between 
common  engineering  measures  for  LO  single-sideband  phase 
noise  and  the  power  spectral  density  of  the  effective  dephasing 
noise  field.  We  unify  familiar  concepts  from  quantum  infor¬ 
mation,  frequency  metrology,  and  precision  oscillator  charac¬ 
terization  in  order  to  translate  phase  noise  specifications  into 
gate  error  probabilities  for  a  variety  of  canonical  single-qubit 
operations  via  the  filter  transfer  function  framework.  Our 
work  highlights  the  fact  that  effective  dephasing  errors  aris¬ 
ing  from  the  phase  fluctuations  in  commonly  used  commercial 
LOs  are  becoming  a  significant  consideration  in  the  context 
of  quantum  information  across  qubit  modalities.  Specifically, 


our  calculations  demonstrate  that  far-from-carrier  phase  noise 
in  state-of-the-art  LOs  provides  a  challenging  upper  bound 
on  operational  fidelities  and  that  the  bandwidth  of  transmit¬ 
ted  thermal  noise  determines  a  non-negligible  contribution  to 
this  error  rate.  Furthermore,  dephasing  errors  due  to  the  high- 
frequency  portion  of  the  phase  fluctuation  spectral  density  is 
not  efficiently  addressed  through  dynamic  error  suppression, 
mandating  improvements  in  LO  hardware  for  long-term  per¬ 
formance  gains. 


II.  RESULTS 

In  a  semiclassical  picture,  qubit  phase  coherence  corre¬ 
sponds  to  maintaining  the  dynamical  phase  of  the  qubit  rel¬ 
ative  to  a  reference  (master)  clock  for  the  system.  Implicit  in 
the  Bloch  sphere  representation  of  a  qubit’s  state  is,  in  general, 
a  transformation  to  a  frame  co-rotating  with  the  nominal  qubit 
Larmor  frequency.  Any  frequency  difference  A  (t)  between 
these  two  effective  oscillators  (the  clock  and  the  qubit)  will 
result  in  phase  accumulation  4>{t)  =  fj  A  (t)dt  in  the  rotat¬ 
ing  frame.  In  a  dephasing  process,  A(t)  evolves  stochastically 
in  time,  introducing  a  degree  of  randomness  to  the  phase  evo¬ 
lution  over  a  time  period  r.  This  simple  observation  indicates 
that  not  only  must  an  experimentalist  consider  the  stability  of 
the  qubit  itself,  but  also  the  stability  of  the  rotating  frame  real¬ 
ized  through  use  of  a  local  oscillator  (LO)  in  the  experiment. 

We  demonstrate  this  observation  formally  by  considering 
an  ensemble  of  identically  prepared  noninteracting  qubits  sub¬ 
jected  to  a  noisy  environment  and  an  external  controller.  Fol¬ 
lowing  our  previous  work  these  dynamics  may  be  incorpo¬ 
rated  by  specifying  the  generalized  time-dependent  Hamilto- 
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FIG.  1.  Schematic  representation  of  the  correspondence  between 
phase  fluctuations  and  extrinsic  dephasing.  A  stable  qubit  interact¬ 
ing  with  a  LO  experiencing  frequency  instability  (captured  by  the 
phase  instablity,  S^lo),  that  is,  the  power  spectral  density  of  phase 
fluctuations),  induces  random  phase  accumulation  in  a  freely  evolv¬ 
ing  qubit,  represented  in  the  Bloch  sphere  picture.  The  same  phe¬ 
nomenology  arises  for  a  LO  outputting  a  perfectly  stable  sinusoidal 
signal  in  the  presence  of  extrinsic  Hamiltonian  terms  that  produce 
uncertainty  in  the  frequency  of  the  qubit  transition. 


nian  as  the  sum  of  control  and  noisy  components 

H(t)  =  Hc(t)  +  H0(t)  (1) 

=  h(f)cr  +  r](t)cr.  (2) 

Here  er  denotes  a  column  vector  of  Pauli  matrices  and  the 
row  vectors  h(i),  77(f)  £  R3  denote  respectively  the  Carte¬ 
sian  components  of  the  control  Hamiltonian,  Hc(t),  and  noise 
Hamiltonian,  H0(t),  in  the  basis  of  Pauli  operators  [12,  13]. 

This  formulation  assumes  an  interaction  picture  co-rotating 
with  the  (resonant)  driving  field.  In  this  picture  the  com¬ 
ponents  h.flt),  i  £  {x,y,z},  of  the  control  field  relate  di¬ 
rectly  to  the  instantaneous  Rabi  rate  f l(t)  =  2||h(i)||,  while 
the  stochastic  noise  fields  rjflt),  i  £  {x,y,  z}  model  semi- 
classical  time-dependent  perturbations  in  each  of  the  three 
spatial  directions.  In  particular,  dephasing  processes  are  cap¬ 
tured  through  the  noise  field  rjz  (t)  describing  stochastic  evo¬ 
lution  along  az . 

In  this  work  we  show  the  formal  equivalence,  in  the  ap¬ 
propriate  interaction  picture,  of  extrinsic  dephasing  noise  pro¬ 
cesses,  and  intrinsic  phase  fluctuations  in  the  control  device 
itself.  That  is,  setting  tjT_v  =  0  to  restrict  attention  to  the  de¬ 
phasing  noise  quadrature,  we  show  the  total  Hamiltonian  takes 
the  form 

H(t)  =  h(t)<r  +  r)z(t)az,  yz(t)  =y(z\t) +rft°\t)  (3) 

where  the  effective  dephasing  term  consists  of  an  extrinsic 
dephasing  field,  and  an  intrinsic  phase-fluctuation- 

induced  field,  on  formally  equal  footing. 

We  model  a  qubit  with  nominal  transition  (angular)  fre¬ 
quency  u>a,  driven  by  a  local  oscillator  with  carrier  frequency 


U4xn  and  write  the  system  Hamiltonian  (h  =  1)  in  the  labora¬ 
tory  frame 


Hs  =-jUacrz  +  -5ua{t)az 

+  0(f)  cos  (ujLOt  +  [ 4>c(t )  +  4>N{t)]jax.  (4) 

The  first  term  corresponds  to  the  Hamiltonian  of  the  qubit 
under  free  evolution,  while  the  second  term  captures  envi¬ 
ronmental  noise  (e.g.  magnetic  field  fluctuations)  effectively 
transforming  the  nominal  energy  splitting  as 

A  +  5ua(fl).  (5) 


The  third  term  corresponds  to  the  qubit-field  interaction  with 
amplitude  Ll(t)  driven  by  a  local  oscillator  with  nominal  car¬ 
rier  frequency  wlo-  and  with  time-dependent  phase  formally 
partitioned  into  a  "control"  component  (desired  phase  evolu¬ 
tion),  4>c{t)i  and  a  “noise”  component  (unwanted  phase  evo¬ 
lution),  </>Ar(f).  Thus,  in  the  lab  frame,  extrinsic  noise  is  cap¬ 
tured  by  5oja(t)  while  intrinsic  noise  in  the  LO  resides  in  the 
LO  phase  fluctuation  field  </>jv(f). 

Moving  to  the  interaction  picture  co-rotating  with  the  car¬ 
rier  frequency  and  making  the  rotating-wave  approximation 
we  obtain 


H 


(<^Lo) 

I 


wLo)dr  +  - 5u>a(t)dz 


<7_ 


} 


(6) 


where  we  define  the  qubit  ladder  operators  a±  =  dx  ±  iay. 

The  LO  phase  fluctuation  produces  an  effective  time- 
dependent  frequency  detuning  through  its  time  derivative, 
SujLO(t)  =  </>Ar(f).  That  is,  the  LO  phase  fluctuation  effec¬ 
tively  transforms  the  LO  frequency  as 


wlo  -a  a>lo  +  8u>ix>(t). 


(7) 


In  order  to  make  this  phenomenology  explicitly  comparable 
with  extrinsic  dephasing,  we  perform  a  second  interaction- 
picture  transformation 


H 


(“LO  ,4>n)  _  rrf  rr(“  Lo) 


=  U  Hj  ’Ui  -H 


<Pn 


(8) 


where  U^N(t)  =  exp[— az\  is  the  evolution  operator 
under  the  LO-noise-induced  dephasing  Hamiltonian  II r:  _:  = 

\4>N(t)dz.  Setting  the  carrier  detuning  to  zero,  uia  —  wlo  = 
0,  the  transformed  system  Hamiltonian  subject  to  noisy  LO 
phase  fluctuations  thereby  takes  the  form 


H 


(cl>lo,0iv) 


-6uja(t)az  -  ~4>N{t)dz 

-0(f)  j  cos[(j}C(t)]ax  +  sin[^c(f)](jy|.  (9) 


Comparing  Eq.  9  with  Eq.  3  we  see  r)i°\f)  =  \ 5uJa(t )  and 


3 


r] ),LO)(i)  =  —i,4>N(t)  =  —^SoJu:)(t),  while  the  third  term  in 
Eq.  9  corresponds  to  the  control  component  h(f)er.  Arriving 
at  this  expression  reveals  the  phase  fluctuation  <pN(t)  in  the 
LO  produces  net  rotations  about  <jz  through  its  time  deriva¬ 
tive,  <j)N{t)-  This  observation  explicitly  links  LO  phase  fluc¬ 
tuations  to  a  dephasing  Hamiltonian  in  a  manner  that  is  in¬ 
distinguishable  from  the  presence  of  an  extrinsic  bath  Hamil¬ 
tonian  Ho(t)  =  r/z(t)az.  This  correspondence  is  captured 
schematically  in  Fig.  II.  Importantly,  once  the  relative  phase 
relationship  between  the  qubit  and  LO  is  established  via  some 
interaction  (the  global  phase  difference  at  t  =  0  is  irrelevant), 
even  for  times  when  VJt)  — »  0  this  Hamiltonian  can  produce 
effective  dephasing  of  the  qubit. 

The  temporal  nature  of  LO  phase  fluctuations  motivates 
the  use  of  the  experimentally  validated  [14-16]  filter  trans¬ 
fer  function  formalism  for  calculating  the  impact  of  this  ef¬ 
fective  dephasing  field  on  qubit  operational  fidelities  (includ¬ 
ing  the  identity  operator  -  memory).  For  a  system  Hamilto¬ 
nian  in  the  form  of  Eq.  3  the  fidelity  of  an  arbitrary  operation 
on  SU( 2)  is  approximated  in  terms  of  the  overlap  integral  of 
Sz  (w),  the  power  spectral  density  (PSD)  of  the  dephasing  field 
r]z{t),  and  GZ!i(lu ),  for  l  £  {x,y,z},  the  filter  transfer  func¬ 
tions  capturing  the  action  of  the  applied  control  Hamiltonian 
Hc(t)  =  h(t)cr  [13,  17], 

In  the  remainder  of  this  paper  we  ignore  extrinsic  dephasing 
by  setting  r]P  (t)  =  0,  thus  restricting  attention  to  dephasing 
induced  by  the  LO.  In  this  case  Sz(io)  refers  to  the  power 
spectral  density  of  the  dephasing  field  rjz  (t)  =  77 iLO)(f)  = 
—  ^</>Ar(f).  However,  the  noise  performance  of  a  LO  is  typ¬ 
ically  characterized  in  terms  of  the  single-sided  phase  noise , 
£(w)  =  \S^(uj)  expressed  in  dBc/Hz,  where  S^(lo)  is 
the  phase  instability,  defined  as  the  unilateral  PSD  of  the  LO 
phase  fluctuations,  (j>N{t)-  In  this  paper  we  include  super¬ 
scripts  (1)  or  (2)  to  denote  unilateral  (one-sided)  and  bilateral 
(two-sided)  PSDs  respectively.  For  ui  >  0  the  unilateral  PSD 
is  simply  twice  the  value  of  the  bilateral  PSD,  representing 
the  double-sideband  spectral  density  in  the  positive  frequency 
domain  (see  Appendix  C2).  To  understand  how  C{lo)  maps 
to  S'!1)  (w),  we  first  define  the  bilateral  PSD  for  phase  fluctua¬ 
tions  (see  Appendix  C  1) 


Co 

jr 

hi 

(I'WI2), 

(10) 

rT 

^t(^)  = 

J  </>N(t)e  iutdt, 

(11) 

where  (•)  denotes  an  ensemble  average  and  the  PSD  is  defined 
symmetrically  over  lo  £  (—00,00).  The  second  line  defines 
the  truncated  Fourier  transform  of  (f>N(t),  where  [— T,  T]  is  an 
experimental  time  window  over  which  phase  fluctuations  are 
measured.  <I>r(a/)  is  an  estimate  of  the  weight  of  the  Fourier 
component  of  <f>N(t)  at  frequency  lo'  ,  generating  sidebands  on 
the  carrier  at  frequencies  wlo  ±  a/.  Let  the  power  in  the  LO 
carrier  be  written  ^(H/7)2,  where  the  coupling  factor  7  <  1 
converts  from  the  amplitude  of  the  LO,  O/7,  to  the  Rabi  rate, 
fi.  Then  the  power  in  S^  (lo)  =  2 S^  (w)  integrated  over  a  1 
Hz  band  centred  on  u> '  is  interpreted  as  the  ratio  of  the  power 


in  these  sidebands  relative  to  4(fi/7)2  (see  Appendix  C4). 

The  phase  noise  C(lo)  is  typically  represented  in  logarith¬ 
mic  units  of  dBc/Hz  as  a  function  of  frequency.  We  denote 
this  by  C(lo)  =  10 log10[£(w)]  to  indicate  logarithmic  units, 
see  Fig.  Ha).  Following  the  above  discussion,  the  Fourier  fre¬ 
quency  lo  in  the  phase  noise  is  offset  from  the  carrier,  with  the 
units  dBc/Hz  indicating  the  sideband  power  value  relative  to 
the  carrier  measured  over  a  1  Hz  band. 

We  convert  from  C(lo)  into  linear  units  via  the  relation¬ 
ship  (oj)  =  2  •  10  to-.  As  demonstrated  above  the 
time-derivative  of  the  phase  is  related  directly  to  a  time- 
dependent  detuning,  and  we  may  relate  their  noise  spectra  as 
=  sP(co)  (see  Appendix  Cl).  The  ultimate 
relationship  therefore  takes  tabulated  single-sideband  phase- 
noise  specifications  to  unilateral  dephasing  power  spectral 
densities  as 


sp  M  =  iW210^.  (12) 

This  correspondence  captures  well-known  [18]  relationships 
between  frequency  and  phase  fluctuations  in  precision  oscilla¬ 
tor  characterization  and  is  represented  graphically  in  Fig.  Ha. 

These  quantities  may  be  directly  incorporated  into  expres¬ 
sions  for  the  fidelity  of  arbitrary  single-qubit  operations  in  the 
presence  of  pure-dephasing  noise  using  [13] 

Fav{r)  «  ^{1  +exp[-x(r)]}, 

V  7  U  is  x,y,z 

=  (i)  E 

'  7  J0  ZS  x,y,z 

^  c",h' 

V  7  0  1 6  x,y,z 

The  transfer  functions  for  the  control,  Gzj(ui),  may  be  calcu¬ 
lated  analytically  and  have  contributions  along  all  Cartesian 
directions  (indexed  by  /),  as  dephasing  noise  present  during  a 
non-commuting  control  operation  (e.g.  Hc  oc  dx)  induces 
both  dephasing  and  amplitude-damping  rotations  [13,  17]. 
The  latter  two  expression  for  x(T)  simply  re-express  the  over¬ 
lap  integral  in  terms  of  various  quantities  encountered  in  de¬ 
scribing  phase  noise. 

We  now  use  these  relationships  to  perform  a  quantitative 
exploration  of  the  impact  of  LO  phase  noise  on  the  fidelity 
of  qubit  operations.  To  calculate  the  operational  infidelity, 
1  —  Favir),  we  insert  the  published  phase  noise  specifica¬ 
tions  for  common  synthesizers  in  the  microwave  X-band,  to¬ 
gether  with  the  analytic  filter  transfer  functions  as  inputs  for 
Eq.  16.  We  consider  two  simple  classes  of  operation:  the 
Identity,  I,  and  a  driven  rotation  by  angle  7r  about  the  a>axis, 
X.  In  each  case  we  calculate  the  relevant  terms  GZii(ut)  fol¬ 
lowing  techniques  in  [13,  16,  17,  19-22]  (see  Appendix  B). 
These  choices  are  representative  of  the  broad  range  of  oper¬ 
ations  required  in  e.g.  the  single-qubit  subset  of  a  universal 


(13) 

(14) 

(15) 

(16) 
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FIG.  2.  a)  Local  oscillator  phase  noise  expressed  as  £(<u)  (left 
axis,  dark  lines)  and  converted  to  S,|1^(w)  (right  axis,  light  lines)  for 
two  grades  of  synthesizer.  Data  extracted  from  specification  sheets 
and  manufacturer  measurements.  For  “lab-grade”  synthesizer  phase 
noise  for  offsets  oj/2n  <  10  Hz  data  are  extrapolated  by  the  man¬ 
ufacturer.  b-c)  Calculated  infidelity  for  different  classes  of  quantum 
logic  operations  (see  main  text).  We  have  numerically  confirmed 
that  for  evolution  times  r  <  100  ms  the  presence  of  a  sharp  low- 
frequency  cutoff  below  uj/2-k  <  1  Hz  has  a  negligible  effect  on  cal¬ 
culations  relative  to  extrapolating  existing  scaling  in  C(u>). 


gate  set.  The  former  allows  insight  into  the  coherent  lifetime 
(often  T2)  of  the  qubit  under  free  evolution  while  the  latter 
informs  how  non-commuting  gate  fidelity  is  reduced  by  the 
presence  of  LO-induced  dephasing  noise. 

We  also  compare  the  performance  of  dynamic  error  sup¬ 
pression  (DES)  strategies  [23-38]  in  reducing  errors  due  to 
typical  phase  noise,  through  application  of  time-dependent 
modulation  of  the  system’s  dynamics  with  the  aim  of  co¬ 
herently  averaging  out  slow  fluctuations.  While  these  pro¬ 
tocols  have  typically  been  associated  with  the  mitigation  of 
environmental  decoherence,  a  growing  body  of  literature  has 
shown  that  errors  induced  by  imperfect  control  -  including 
LO  phase  noise  -  can  also  be  mitigated  using  the  same  tech¬ 
niques  [4,  14,  16].  In  our  case,  the  analytic  form  for  the  DES 
version  of  the  identity  is  given  for  spin  echo,  while  for  a  driven 
rotation  by  angle  7r  (again  ex  ax),  we  employ  a  form  defined 
as  both  a  dynamically  corrected  gate  [34—36]  and  a  Walsh  am¬ 
plitude  modulated  filter  [16,  39].  Full  analytic  expressions  for 
these  filter  functions  may  be  found  in  Appendix  B. 

For  this  analysis  we  have  selected  two  distinct  grades  of  LO 
to  demonstrate  the  significance  of  LO  choice  in  quantum  con¬ 
trol  experiments.  Figure  Ha  shows  the  single-sideband  phase 


noise  at  10  GHz  for  both  a  “lab  grade”  synthesizer  (Vaunix 
LMS-123)  and  a  “precision”  synthesizer  (Keysight/ Agilent 
8267D  OPT-UNY).  Across  much  of  the  band,  as  expected, 
the  phase  noise  for  the  precision  synthesizer  is  approximately 
40  dB  lower  than  the  lab-grade  unit.  For  the  synthesizers  stud¬ 
ied  here  C(ui)  declines  rapidly  with  offset  frequency  from  the 
carrier,  exhibiting  various  power-law-dependences.  Details  of 
the  origins  of  the  underlying  processes  may  be  found  in  [40]. 

Calculations  of  infidelity  for  the  four  operations  outlined 
above  are  presented  in  Fig.  Ilb-c  and  summarized  in  Table  I. 
For  the  identity  operator  we  find  an  approximate  improvement 
of  104  in  residual  error  rate  due  to  use  of  a  high-precision 
frequency  source,  and  the  infidelity  remains  below  10  "  4  out 
to  100  ms  evolution  time.  By  contrast  the  lab-grade  syn¬ 
thesizer  induces  an  error  exceeding  ~  0.1%  beyond  a  few 
microseconds  of  evolution  time.  These  findings  are  relevant 
for  trapped-ion  [41],  superconducting  [42],  and  semiconduc¬ 
tor  [43]  qubits.  The  plateau-like  behavior  in  infidelity  with 
increasing  evolution  time  is  due  to  the  interplay  of  the  power 
spectrum  and  filter-function  [44],  and  is  similar  to  phase- 
error  saturation  phenomenology  observed  in  precision  oscil¬ 
lator  characterization  [40]. 

In  both  free  evolution  and  driven  operations  the  deleterious 
impact  of  using  a  lab-grade  synthesizer  increases  as  the  du¬ 
ration  of  the  operation  grows.  For  instance,  for  driven  opera¬ 
tions  as  long  as  100  /is,  particularly  relevant  for  atomic  qubits, 
the  precision  synthesizer  only  induces  an  error  ~  1 0  "  ‘  while 
the  lab-grade  synthesizer  induces  an  error  more  than  20, 000  x 
larger,  reaching  ~  0.2%.  The  primary  reason  for  this  behav¬ 
ior  is  the  enhanced  far-from-carrier  phase-noise  performance 
in  the  precision  synthesizer,  especially  over  the  band  from  ap¬ 
proximately  1  —  1000  kHz  offset. 

An  interesting  observation  is  that  in  the  presence  of  these 
noise  power  spectral  densities  and  the  realistic  evolution  times 
selected  for  our  calculations,  DES  protocols  have  a  minimal 
effect  on  gate  performance.  Only  for  very  short  driven  op¬ 
erations  (<100  ns)  does  DES  provide  performance  enhance¬ 
ments,  where  for  gates  with  r  <  20  ns  error  rates  are  sup¬ 
pressed  more  than  1000  x.  By  contrast,  implementing  the 
identity  operator  over  most  of  the  range  r  £  {1/is,  100ms} 
the  use  of  spin  echo  actually  degrades  operational  fidelity. 

This  is  explained  by  considering  that  while  C(oj)  appears 
to  decline  rapidly  with  offset  frequency,  transforming  this 
spectrum  to  (Fig.  Ha,  right  axis)  reveals  the  high- 

frequency  dominance  of  the  resulting  dephasing  noise,  reflect¬ 
ing  an  approximately  Olimic  spectrum  (Sz  oc  w)  [45,  46],  In 
this  regime  DES  is  known  to  perform  inefficiently  as  the  noise 
evolves  rapidly  compared  to  the  control  and  the  physics  of  co¬ 
herent  averaging  fails:  a  violation  of  the  so-called  decoupling 
limit. 

Thermal  noise  in  oscillators  imposes  a  lower-error  bound  to 
be  considered;  even  if  LO  hardware  were  improved  we  could 
do  no  better  than  saturating  the  thermal  noise  floor  across  the 
control  bandwidth.  This  is  generically  quoted  as  -174  dBm/Hz 
for  a  matched  load  at  290K.  This  is  an  absolute  noise  power, 
meaning  that  for  LOs  with  power  0  dBm  the  single-sideband 
phase  noise  floor  will  take  a  value  Cmin  =  —174  dBc/Hz, 
and  it  will  rise  (worsen)  with  decreasing  LO  power.  Higher- 
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FIG.  3.  Calculated  infidelity  floor  (fidelity  ceiling)  imposed  by  ther¬ 
mal  noise  at  4K  and  290K  with  variable  control  bandwidths.  This 
calculation  is  independent  of  operation  length,  r  under  the  condi¬ 
tions  outlined  in  the  main  text.  Results  hold  for  Ramsey  and  both 
primitive  and  WAMF  7r  rotations.  Infidelity  floor  for  spin  echo  is 
three  times  larger.  See  Appendix  D 

power  LOs  can  in  principle  have  lower  phase-noise  floors, 
but  will  generally  suffer  from  additional  amplifier-induced 
phase  noise.  We  therefore  consider  the  thermal  noise  floor 
for  a  LO  at  0  dBm  output  as  an  approximate  lower-bound 
for  the  phase  noise  and  calculate  the  resulting  upper-bound 
on  operational  fidelity  in  the  presence  of  this  noise.  As¬ 
suming  this  constant  minimal  value  for  the  phase  noise  over 
frequency,  C(u>)  =  Cmin,  the  infidelity  is  computed  using 
Xmin  ~  (kwc/ 2^^)10Cr,lirll,10  (see  Appendix  D).  In  this  ex¬ 
pression  u)c  is  the  cutoff  frequency  for  the  control  bandwidth, 
k  is  a  characteristic  scaling  factor  depending  on  the  control 
protocol,  and  Xmin  is  approximately  independent  of  r  in  the 
limit  ujc/2i r  t~1. 

This  straightforward  calculation  reveals  that  broadband 
noise  at  room  temperature  imposes  a  non-negligible  upper 
bound  on  gate  fidelity  as  shown  in  Fig.  3.  For  a  typical  20  GHz 
bandwidth  associated  with  standard  coaxial  cable  connections 
between  LO  and  qubit,  the  thermal  noise  floor  induces  qubit 
infidelity  in  excess  of  10-8.  Reducing  the  thermal  noise  floor 
to  that  associated  with  a  4K  bath  (-192  dBc/Hz)  improves  the 
fidelity  by  two  orders  of  magnitude.  Similarly,  restricting  the 
control  bandwidth  to  10  MHz,  e.g.  via  narrowband  filtering 
around  the  carrier,  sets  the  infidelity  <  10-11. 


Superconducting 

Trapped  Ion 

Lab-Grade 

Primitive 

4.1  x  10~8 

7.7  x  10"a 

WAMF 

1.5  x  10-ii 

2.3  x  10 b 

Precision 

Primitive 

1.9  x  10-1U 

8.3  x  10"8 

WAMF 

3.4  x  10-i8 

4.3  x  10"s 

290K  (4K)  Floor 

100  MHz 

6.3  x  10 11(1  x  10~i:i) 

1  GHz 

6.4  x  10~10  (1  x  10"11) 

10  GHz 

6.4  x  10-y  (1  x  10-iu) 

TABLE  I.  Comparison  of  typical  error  rates  for  superconducting 
(10  ns)  vs  trapped-ion  (10  p, s)  driven-gate  operations.  Main  re¬ 
sults  calculated  assuming  an  integration  bandwidth  of  10  MHz.  The 
temperature-dependent  noise-floor  for  different  cutoff  frequencies  is 
also  shown  -  these  data  are  largely  independent  of  operation  time  and 
type.  Noise  floor  calculations  strictly  valid  only  for  cutoff  frequen¬ 
cies  high  compared  with  the  inverse  operation  time. 


m.  DISCUSSION 

The  calculations  we  have  presented  demonstrate  that  far- 
from-carrier  phase  noise  dominates  residual  errors  for  com¬ 
mercial  sources  and  is  not  easily  suppressed  by  DES  pro¬ 
tocols  due  to  the  high-frequency  weight  in  the  dephasing 
power  spectrum.  Our  results  enhance  existing  arguments 
about  the  merits  of  qubit  modalities  that  accommodate  fast 
control  pulses;  for  instance,  reducing  driven  operation  times 
from  10  p s  to  10  ns,  upper-bound  fidelities  can  be  improved 
by  more  than  five  orders  of  magnitude  for  a  high-performance 
synthesizer.  Interestingly  these  example  figures  represent 
common  operation  times  for  atomic  vs  solid-state  qubits  re¬ 
spectively. 

However,  any  move  to  shorter  control  operations  must  be 
balanced  with  cognizance  of  the  control  bandwidth.  For  ex¬ 
ample,  we  indicated  that  bandwidth-dependent  thermal  noise 
floors  may  pose  appreciable  upper  bounds  on  operational  fi¬ 
delities.  Given  the  exponential  scaling  of  infidelity  with  con¬ 
trol  bandwidth,  there  is  a  tradeoff  between  modulation  band¬ 
width  and  achievable  error  rates  that  must  be  considered  - 
solely  maximizing  pulse  rates  may  not  ultimately  yield  op¬ 
timal  performance.  In  addition,  the  impact  of  suppressing  the 
thermal  noise  floor  from  -174  dBc/Hz  to  -192  dBc/Hz  sug¬ 
gests  that  cryogenic  LOs  [47]  embedded  in  cryogenic  control 
architectures  [48]  and/or  careful  cryogenic  engineering  of  am¬ 
plified  room-temperature  LOs  will  play  an  important  role  in 
engineering  high-fidelity  quantum  information  systems.  We 
believe  it  is  likely  these  considerations  will  have  significant 
bearing  on  how  qubits  are  connected  to  classical  control  elec¬ 
tronics  and  that  a  detailed  cost-benefit  analysis  will  need  to  be 
performed  weighing  operation  and  modulation  rates  against 
added  broadband  phase  noise. 

Our  observations  are  particularly  concerning  because  they 
imply  that  in  many  cases  clock-induced  errors  will  pose  a 
lower  bound  on  error  rates  feeding  into  quantum-error  cor¬ 
rection  and  fault-tolerance  calculations  without  hardware  im¬ 
provement  in  a  technology  space  which  has  received  some¬ 
what  limited  attention.  Using  precision  LOs  in  most  ex¬ 
perimental  labs  will  suppress  otherwise  performance-limiting 
error  rates  below  currently  measurable  limits.  However, 
the  bounds  on  error  rates  associated  with  both  free  and 
driven  evolution  are,  in  the  view  of  these  authors,  remark¬ 
ably  high.  Given  the  efficacy  and  resource-efficiency  of  QEC 
is  sensitively  dependent  to  how  far  the  error  rate  is  below 
the  fault-tolerance  threshold  [49],  error  floors  ~  10-'  for 
microsecond-scale  evolution  times  using  precision  LOs  make 
clear  the  need  to  focus  on  developments  in  advanced  master- 
clock  hardware  in  addition  to  qubit  improvements. 

While  we  have  focused  on  phase  fluctuations  on  short 
timescales,  consideration  of  both  fast  phase  fluctuations  and 
long-term  LO  stability  is  required.  It  is  important  to  distin¬ 
guish  the  impact  of  long-term  instabilities  (due  to  close-to- 
carrier  phase  fluctuations)  on  experimental  measurements  of 
qubit  fidelities  from  “intrinsic”  error  rates;  the  process  of  data 
acquisition  and  averaging  over  many  individual  experiments 
contributes  error  due  to  slow  drifts  that  dominate  the  actual 
error  rate  experienced  in  any  individual  operation  (see  [22] 
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for  details).  In  addition,  long-time  instabilities  may  be  detri¬ 
mental  to  systems  with  multiple  control  generators,  particu¬ 
larly  if  the  generators  exhibit  phase  diffusion  with  respect  to 
one-another.  A  detailed  accounting  of  such  architectural  is¬ 
sues  will  be  the  subject  of  future  work. 

In  summary  our  results  have  unified  concepts  from  fre¬ 
quency  metrology  and  quantum  information,  permitting  ex¬ 
perimentalists  to  directly  translate  common  LO  hardware 
specifications  into  meaningful  estimates  of  qubit  coherence 
and  operational  fidelities.  We  have  discovered  that  far-from- 
carrier  phase  noise  and  broadband  thermal  noise  pose  a  per¬ 
formance  limiting  upper-bound  on  operational  fidelities,  and 
improving  these  bounds  will  require  major  hardware  improve¬ 
ments.  As  a  result  we  expect  to  see  a  growing  emphasis  on 
high-performance  LO  synthesis  chains  in  experimental  quan¬ 
tum  information  laboratories  across  all  technology  platforms, 
and  growing  importance  of  LO-induced  errors  in  quantum 
systems  design.  Future  studies  will  examine  not  only  clock 
synthesis,  but  clock  distribution,  with  an  eye  towards  archi¬ 
tectural  impacts  of  clock  distribution  in  quantum  systems. 
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Appendix  A:  Phase  noise  in  LO  upconversion 

An  experimental  LO,  operated  near  resonance  with  the 
qubit  Larmor  frequency  will  generally  be  upconverted  from 
a  lower-frequency  phase  reference.  This  process  results  in  ad¬ 
ditional  phase  noise;  assuming  noiseless  upconversion  by  a 
factor  of  N  to  reach  the  desired  carrier,  the  carrier  phase  noise 
in  log  units  £cairier(w)  -A  £Ref(w)  +  20 log N  dB  [18].  The 
second  term  on  the  right  indicates  that  the  frequency  upcon¬ 
version  leads  to  an  effective  increase  in  the  fractional  phase 
instability  of  the  carrier.  The  synthesis  chain  will  also  gener¬ 
ally  add  some  phase  noise  on  top  of  this  quantity. 


Appendix  B:  Analytic  forms  of  filter  functions 


function  component  takes  the  form 


Here  we  give  the  explicit  form  for  the  filter-transfer  functions 
for  the  qubit  operations  detailed  in  the  main  text.  Namely, 
Ramsey  (free  evolution),  spin  echo,  primitive  7r-pulse,  and 
DCG/WAMF  7r-pulse.  In  the  dephasing  noise  quadrature,  the 
filter  function  in  all  cases  takes  the  generic  form 

Fz{  w)=  GM  (Bl) 

x,y,z 

where  the  transfer  functions  for  the  control,  Gz  i(u>),  capture 
contributions  along  all  Cartesian  directions  (indexed  by  l ).  In 
general  all  three  components  are  non-zero  as  dephasing  noise 
induces  both  dephasing  and  amplitude -damping  rotations  [13, 
17]  during  non-commuting  control  operation  (e.g.  Hc  cx  ax). 

In  the  following  sections  we  give  the  functional  form  only 
for  the  nonzero  transfer  functions  G  zj  (lu  )  for  each  of  the  qubit 
operations  listed  above. 

1.  Ramsey  (Spin  Echo) 

For  free  evolution  (i.e.  zero  driving  field)  over  a  duration  r 
the  only  nonzero  transfer  function  component  takes  the  form 
[15,  19-21] 

G£an-ey)  (w)=4sin2(0).  (B2) 

2.  Spin  Echo 

Spin  echo  is  the  DES  version  of  the  identity  operation.  This 
takes  the  form  of  free  evolution  over  duration  r,  but  with  an 
instantaneous  nx  rotation  at  time  r/2.  The  nonzero  transfer 


Ggcho)(w)  =  1  +  eiuT  -  2eiwT/2  2  =  16  sin4  (0) . 

(B3) 


3.  Primitive  7r-pulse 

For  a  finite-strength  rotation  of  duration  r  and  Rabi  rate 
O  =  7 r/r,  there  are  two  nonzero  transfer  functions: 

.  .2  2 

G&-prim)M=  ^3^(eiWT  +  1)  -  <B4) 

GirPnm)M  =  |o^t(e^  +  l)|2.  (B5) 

The  second  line  arises  because  of  the  noncommuting  nature 
of  the  control  (cx  ax )  and  the  dephasing  noise  (cx  az),  as  de¬ 
scribed  in  [13,  17]. 

4.  DCG/WAMF  7r-pulse 

This  is  a  finite-strength  sequence  performing  a  net  n  rotation 
and  defined  both  as  a  dynamically  corrected  gate  (DCG)  and 
a  Walsh  amplitude  modulated  filter  (WAMF)  [16,  39],  Let 
Rx(0 ,r)  be  a  non-instantaneous  qubit  rotation  about  ux  on 
the  Bloch  sphere,  through  an  angle  9  over  a  duration  r  (i.e. 
with  Rabi  rate  O  =  0/t).  Then  the  DCG/WAMF  sequence  of 
total  duration  r  is  defined  by  the  operator  sequence 

Rx  (tt,  0  Rx  0,  0  Rx  (n,  0  ,  (B6) 

with  transfer  function  components 
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where  \£\  is  the  number  of  elements  in  the  set  8.  We  include 
explicitly  the  subscript  k  in  the  expression  for  the  time  average 
to  emphasize  that  it  is  computed  over  a  single  realization  of 
y(t).  We  assume  the  following  properties  of  y(t): 

1.  y(t)  is  a  wide-sense  stationary  (w.s.s.)  process. 

2.  y(t)  is  an  ergodic  process. 


3.  y(t)  is  a  zero-mean  process:  (y(t))  =  yk(t)  =  0,  Vi. 

By  w.s.s.  we  mean  the  ensemble  autocorrelation  function 
Cy(r)  is  invariant  under  time-translations  of  r.  By  ergodic 
we  mean  that,  in  the  limit  of  infinite  ensemble  sizes,  the  en¬ 
semble  mean  of  y(t)  -  or  the  ensemble  mean  of  a  function  of 
y(t)  -  approaches  the  corresponding  time-average  over  any 
given  realization  yk(t).  That  is, 

(. f{y{t )))  * — >  f{Vk{t))-  (C3) 

In  particular,  ergodicity  permits  the  autocorrelation  function 
to  be  equivalently  defined  either  in  terms  of  the  ensemble  av¬ 
erage  or  the  time  average 

0/(r)  =  (y(t)y(t  +  t))  =  yk{t)yk{t  +  r).  (C4) 

Due  to  the  assumption  of  w.s.s.  the  autocorrelation  function 
above  depends  only  on  the  time  difference,  r,  between  mea¬ 
surements.  We  now  define  the  truncated  Fourier  transform  of 

y(t)  by 


YT(to)  =  j  ^y(t)e~iutdt,  (C5) 

where  the  angular  frequency,  non-unitary  convention  for 
Fourier  transforms  has  been  used.  This  definition  avoids  the 
difficulty  that  the  full  Fourier  transform  may  not  be  well  de¬ 
fined  if  the  integral  over  y(t)  fails  to  converge  in  the  limit 
T  — >  oo.  The  power  spectral  density  (PSD)  of  y(t)  may  then 
be  defined  by 


Sy{u)  =  ^  (lyr(w)|2)  •  (C6) 

Writing  this  out  we  explicitly  we  obtain 

Sv(u)  =  ^ 

(Cl) 

=  lim  ±  [T  [T  (ymt^e^-^dtdf 

T-> oo  2 T  J_T  J_T 

(C8) 

=  lim  f  [  Cy(r)e~iuT dtdt' ,  (C9) 

T— ^oo  11  J _T  J _T 

where  we  have  used  the  fact  that  y(t)  is  real,  and  defined  the 
new  variable  t  =  t'  —  t.  Transforming  the  integral  domain  we 


get 


=  lim  -J-  [ ^  (2 T  -  \r\)Cy(r)e-^dT  (CIO) 

=  lim  /  Cy(T)e~iuTdT  (Cll) 

T->oo  J _2T 

+  lim  tL  [  M  Cy(T)e~lWT dr.  (C12) 

Assuming  y(t)  has  finite  temporal  correlations,  appropriate 
for  realistic  physical  processes,  the  autocorrelation  function 
Cv(t)  vanishes  in  the  limit  r  — >  oo.  Consequently  the  second 
term  above  vanishes  in  the  limit  T  — ►  oo  yielding 

/OO 

Cy(T)e~iuiTdT ,  (C13) 

-oo 

Cy(T)  =  ^~  [  Sy(u)eiuiTdi0.  (C14) 

2n  J- oo 

That  is,  Sy  (ui)  and  Cy  (r)  form  a  Fourier  transform  pair,  which 
is  a  statement  of  the  Wiener-Khintchine  theorem. 


2.  Bilateral  and  Unilateral  PSDs 

The  power  spectral  density  formulated  in  Eq.  C6  is  by  defi¬ 
nition  a  non-negative  even  function  on  the  real  domain  ut  £ 
(—oo,  oo).  This  is  known  as  the  bilateral  PSD,  with  total 
power  (or  mean-square  fluctuation)  given  by  integrating  over 
both  positive  and  negative  frequency  domain: 


A*  =  0,(0)  =  < y(t )2> 

(C15) 

=  2 \  j 

(C16) 

=  /  Sv[y)dv ,  u>  =  2itv. 

J  —oo 

(C17) 

The  bilateral  PSD  naturally  arises  in  mathematical  analy¬ 
sis  involving  Fourier  transformations,  satisfying  the  Wiener- 
Khinchine  theorem  describing  its  relationship  with  the  auto¬ 
correlation  function  Cv(t)  in  Eqs.  C13  and  C14. 

It  is  also  useful  when  considering  systems  subject  to  quan¬ 
tum  noise,  for  which  the  corresponding  quantum  operator  y(t) 
in  general  does  not  commute  at  different  times  and  is  in  princi¬ 
ple  a  complex  number,  i.e.,  (y(t)y(t  +  r))  =  a  +  ib.  The  spec¬ 
tral  density  of  a  complex  number  is  no  longer  an  even  function 
in  !jj  so  must  be  specified  for  both  positive  and  negative  fre¬ 
quencies.  Such  non-symmetric  spectra  arise,  for  example,  in 
the  context  of  qubit  stimulated  and  spontaneous  emission  and 
absorption  to/from  the  environment. 

For  real-valued  y(t),  however,  the  bilateral  PSD  is  sym¬ 
metric  and  may  be  replaced  by  the  unilateral  PSD,  more  typ¬ 
ically  used  in  engineering  and  metrological  applications  deal¬ 
ing  with  manifestly  classical  signals.  Environmental  fluctu¬ 
ations  responsible  for  qubit  dephasing  fall  into  this  category. 
For  clarity,  in  this  paper  we  use  the  following  notation  to  keep 
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the  usage  of  bilateral  (two-sided)  and  unilateral  (one-sided) 
PSDs  distinct: 

S^2\uj)  ( bilateral ),  (C18) 

( unilateral ),  (C19) 


Using  the  definition  of  <f>r(w),  the  fact  that  5u>ix){t)  =  <pN{t), 
and  integrating  by  parts  we  find 

WT(w)  =  +  R,  (C32) 

R  =  <fiN(T)e~iuiT  -  ^iv(-T)eiwT,  (C33) 


where  the  unilateral  PSD  is  defined  by 


S{1)  AA  =  f  242)  M  for  w  >  0 
y  [  ’  [0  for  w  <  0 


(C20) 


With  these  definitions  the  total  power  in  the  signal  is  given  by 


Ptot  — 


(C21) 

(C22) 

(C23) 


where  the  residual  term  R  comes  from  the  integration  limits 
after  integrating  by  parts,  and  we  have  used  the  Fourier  re¬ 
lation  d/dt  iui.  Taking  the  modulus  square  we  therefore 
obtain 


|WV(w)|2  =w2|$tM|2  +  P',  (C34) 

where  the  residual  terms  are  absorbed  into 

R'  =  \R\2  +  [)w$T(w)f?*]  +  Hw$T(w)*f?] .  (C35) 


Consequently 


SWo M  =  Tlim_  ^  (u2  |$t(w)|2  +  R'^  (C36) 

rlim  ^(|$t(w)|2)  +  lim  (R!) . 

V— yoc  Z1  \  /  T— yea  Z1 


3.  PSDs  for  Phase  and  Frequency  Noise 

Consider  a  local  oscillator 

f 1(f)  cos[wLof  +  4>N{t)\  (C24) 

with  instantaneous  frequency 

w(t)  =  —  [tULof  +  <^jv(f)]  =  WLo  +  <^Ar(f),  (C25) 

or 

w(f)  =  wlo  +  ^WLo(f)  (C26) 

where  the  (angular)  frequency  noise  takes  the  form 

Suva®  =j>N{t).  (C27) 

Following  Eq.  C5  we  define  the  truncated  Fourier  transforms 
for  phase  and  frequency  noise  fields 

$T(w)  =  J  <j)N{t)e-iutdt,  (C28) 

Wt(w)  =  J  SujLO{t)e-iutdt.  (C29) 

Then  the  phase  and  frequency  noise  PSDs  are  defined  by 

SWM  =  T1™0^(l$T(a;)l2)’  (C3°) 

55wlo(w)  =  ^lim^  ^(\WT(u)\2y 


The  term  in  square  brackets  is  the  definition  of  Scf,N(uj).  The 
second  term,  on  the  other  hand,  vanishes  in  the  limit  T  — >  oo 
as  the  residual  term  R'  does  not  grow  sufficiently  rapidly  with 
T.  We  therefore  set 

lim  <#>  =  0  (C38) 

T  — ^oo  A jL 

from  which  it  follows 

SSuu>(u)  =w2S0JVH.  (C39) 

Now  define  the  noise  field 

/3(f)  =  adujLO(t)  =  a</>N(t)  (C40) 

by  simply  scaling  the  frequency  noise  by  some  factor  a.  Sub¬ 
stituting  this  into  the  above  derivation  we  obtain  the  slightly 
more  general  PSD  relation 

S0(u)  =  a2SsUw  =  a2  u2S (C41) 

These  relations  apply  to  bilateral  and  unilateral  PSDs. 


4.  Interpretation  of  S<j,N  (w) :  Sideband  and  Modulation 
Theory 

Here  we  provide  an  interpretation  of  the  phase  instability 
S<pN  (u>)  defined  in  Eq.  10  as  the  ratio  of  the  power  in  the 
sidebands  relative  to  the  power  in  the  carrier.  Consider  the 
LO  signal 


(C31) 


V(t)  =  Qcos[uix>t  + 


(C42) 
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where  in  a  given  realization,  over  some  window  [—T,T], 
(f>]\f(t)  has  spectral  components  described  by  some  Fourier 
spectrum.  To  obtain  a  simple  interpretation  of  the  resulting 
PSD  S(f,N(uj),  we  restrict  attention  to  the  case  where  4>n(1) 
consists  of  a  single  (phase-randomized)  Fourier  component 

</>Ar(f)  =  +  c.c.)  =  oij  cos (ujjt  +  ipj),  (C43) 

equivalent  to  imposing  a  narrow-band  filter  at  Fourier  fre¬ 
quency  ajj.  The  resulting  waveform  may  be  expressed  in 
terms  of  an  infinite  comb  of  frequency  components, 

OO 

V(t)  =  fl  ^2  Jn(&j)  cos[(cjlo  +  nujj)t  +  nipj],  (C44) 


weighted  by  Bessel  functions  Jn(x)  and  centred  around  the 
carrier  frequency  wlo-  Assuming  the  phase  fluctuations  are 
small  relative  to  the  carrier  frequency  (ay  <C  0J\_o)  we  may 
make  the  small  angle  approximation  to  obtain 

V(t)  «  V0(t)  +  V+(t)  +  (C45) 

where 

Vo(f)  =  flcos(tULof)  ( LO/carrier ),  (C46) 

V+(t)  =  -^oij  cos [(cu+t  +  ipj)]  (right  sideband),  (C47) 

n 

V-(t)  =  —  cty  cos[(tu_f  —  ipj)\  (left  sideband),  (C48) 

and  ia±  =  ojlo  ±  ui7  are  the  positive-  and  negative-detuned 
frequencies  of  the  principle  (first-order)  sidebands  on  the  car¬ 
rier.  Using  Eqs.  C5  and  C6  we  may  now  compute  the  PSDs 
for  Vo(t)  and  V±(t). 


Substituting  Eq.  C43  into  Eq.  C5,  the  truncated  Fourier 
transform  of  the  phase  fluctuation  is  given  by 

$t(w)  =  ctjT  (e'^sinc [T(w  —  u;y)]  +  e_J^3sinc [T(w  +  tuy)])  . 

Using  Eq.  C6  the  (bilateral)  PSD  is  therefore  given  by 

(w)  =  lim  — 

oo  2  TV  / 

To? 

=  lim  -  (sinc2[T(w  —  cuy)]  +  sinc2[T(w  +  wy)])  , 
T—foo  2 

where  cross  products  in  the  modulus  square  between  the  dif¬ 
ferent  frequency  components,  sinc[T(w  ±  ojy )] ,  are  omitted 
as  these  vanish  in  the  limit  T  — >  oo.  Carrying  out  the  limit, 
subject  to  the  condition  that  the  area  under  the  sine2):/:)  curve 
remains  equal  to  n,  we  obtain 

S™  M  =  \ -  wy )  +  S (w  +  wy )] .  (C49) 


Similarly, 

Svo(u)  =  ~  wLo)  +  +  wlo)]  ,  (C50) 

=  \  [5(tu-u;±)-t-5(a;  +  a;±)].  (C51) 


Let  Py(u},  A)  denote  the  power  in  the  signal  y(t)  in  (angu¬ 
lar)  frequency  bandwidth  A  centred  on  (angular)  frequency 
w,  namely  over  bandwidth  [tu  —  A/2,  ut  +  A/2].  Then 


Pv(u),A) 


(C52) 


where  the  factor  of  2  inside  the  integral  accounts  for  the  equal 
contribution  to  the  total  power  in  this  bandwidth  from  the  neg¬ 
ative  frequency  domain  of  the  bilateral  PSD.  In  our  case,  for 
any  A  >  0,  we  obtain 


P&N  (wy  >  A)  =  ^  a2 ,  (C53) 

PVo(coLOl  A)  =  ^U2,  (C54) 

Pv±(u±,  A)  =  ^(^)2-  (C55) 


The  ratio  of  the  total  power  in  both  sidebands  relative  to  the 
carrier  is  therefore  given  by 


Py+  (tu+ ,  A)  +  Pv_  (tu_ ,  A) 
Pvai^LO,  A) 


P*w(wy,A).  (C56) 


In  particular,  setting  A  =  1/27T  rad  =  1  Hz  and  dividing  the 
power  by  the  measurement  bandwidth,  we  obtain 


p<l>N(uj,A)  _  ^2  [rad2 / rad] .  (C57) 


Using  Eq.  C49,  this  compares  directly  with  the  (unilateral) 
PSD  for  4>N(t)  (i-e-  the  phase  instability) 


=  tt a2<5(w  -  uij).  (C58) 


Consequently,  S (ui)  may  be  interpreted,  according  to  Eq. 
C56,  as  the  ratio  of  the  power  in  the  sidebands  relative  to  the 
power  in  the  carrier  measured  over  a  1  Hz  bandwidth. 


Appendix  D:  Thermal  noise  floor 


Here  we  consider  the  contribution  to  infidelity  from  the  funda¬ 
mental  limit  associated  with  a  thermal  noise  floor.  The  theo¬ 
retical  minimum  noise  power  for  an  oscillator  coupled  into  50 
O  cable  is  kTB  for  bandwidth  B,  temperature  T,  and  where 
k  is  the  Boltzmann  constant.  Normalized  to  a  1Hz  bandwidth 
this  is  expressed  in  dBm/Hz  (Table  II). 
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T[K] 

Cmin  [W/HZ] 

Cmin  [dBm/Hz] 

4 

4  x  10"21 

-174 

290 

4  x  10"21 

-192 

TABLE  II.  Lower  bounds  on  phase  noise  set  by  thermal  noise  floor. 


We  assume  an  oscillator  at  0  dBm  where  the  minimal  phase 
noise  is  given  by  the  power  in  the  thermal  noise  referenced  to 
the  power  in  the  carrier.  We  impose  a  high-frequency  cutoff 
loc  on  to  the  otherwise  white  thermal  noise  to  account  for  fil¬ 
tering  due  to  e.g.  coaxial  cable  cutoff  frequencies  or  narrow- 
band  filtering.  Hence,  the  thermal  noise  floor  may  be  charac¬ 
terized  by 


C(u) 


£-"min  0  —  W 
0  LO  >  UIC 


(Dl) 


The  infidelity  1=1  —  T  associated  with  this  noise  floor  is 
given  by 


1.  Ramsey  (Free  Evolution) 


In  the  case  of  free  evolution  we  have 


rUc  /ujt\2 

dujFz(uj)  =  J  4sin^  —  j  duj 
sin  (wct) 


=  2wc  (  1  - 
~  2cuc, 


U)CT 
UJCT  »•  1. 


Consequently 


1  Cmin 

x(r)  =  ^-!0  10  (2  0JC)  (Ramsey) 

Z7T 


(Dll) 

(D12) 

(D13) 


(D14) 


is  proportional  to  the  cutoff,  with  proportionality  factor  as¬ 
sociated  with  strength  of  the  thermal  noise  floor.  Using  the 
appropriate  filter  function,  similar  results  hold  for  other  qubit 
operations. 


1=1-  - 
2 


1  +  e"x(T) 


where 


X(r)  = 


l 

1  £„ 

—  10^ 
27 r 


/o 


(D2) 


)Cw  Fz{u) 

(D3) 

■OJc 

dcoFz  (to) 

(D4) 

2.  Spin  Echo 


For  the  spin  echo  sequence  it  is  just  as  straightforward  to  de¬ 
rive 


and  Fz(io)  =  x  z  Gzj(u>)  is  the  filter  function.  Thus, 
the  infidelity  is  given  by  the  integral  of  the  filter  function  over 
the  band  [0,  wc]  scaled  by  a  constant  factor  associated  with  the 
strength  of  the  thermal  noise  floor.  For  the  qubit  operations 
described  in  the  main  text  we  can  obtain  analytic  expressions 
for  this  integral.  We  assume  the  cutoff  frequency  ojr:  is  fast 
compared  to  the  duration  r  of  the  qubit  operation: 


du>Fz(oj ) 


6wc. 


Hence 


X(r) 


27T 


(6wc) 


( spin  echo). 


(D15) 


(D16) 


1 

0JC  »  -  <*=>  UJCT  »  1. 
T 


(D5) 


In  this  case  the  bulk  of  the  thermal  noise  floor  lies  in  the  band¬ 
width  beyond  the  stopband  of  the  filter  function  and  the  inte¬ 
gral  is  very  well  approximated  by  its  asymptotic  value,  linear 


in  the  cutoff  frequency.  Specifically 

X(r)  =  |JlO%%  (D6) 

where 

k  =  2  (Ramsey),  (D7) 

k,  =  6  (spin  echo),  (D8) 

k  =  2  (primitive  n-pulse),  (D9) 

k  --  2  (DCGfWAMF  it -pulse).  (DIO) 


3.  Primitive  7r-pulse 


The  form  of  the  filter  function  gets  rapidly  more  complicated 
for  finite-width  pulses.  For  the  primitive  7r  pulse  we  have 


2t tV  2  z4 

=  (n*-z*)*  +  (n*-z*)* 
2n2z2cos(z)  2z4  cos(z) 

+  (^-*2)2  +(7r2-^)2> 


(D17) 


where  we  introduce  the  dimensionless  variable  z  =  lot.  The 
main  contribution  to  the  integral  over  the  band  [0,wc]  for 
loc  1/t  comes  from  those  terms  in  Eq.  D17  which  sur- 


We  derive  these  expressions  below. 


vive  the  limit  2  — >  00.  We  find 


r 

hm  -r~2 - 2T2  = 

2-^00  —  zzy 

2z 4 

lim  T~ 9 - 2t2  =  2> 

z-)-oo  (7 Tz  —  Zz)z 

2tt2z2cos(z ) 


z-t 00  (7T2  —  22)2  [z->oo  (7T2  —  22) 

224  cos(2)  I"  2z4 

lim  t ^ — 2W  ^  lim  — 2W 

z-t 00  (7T-  —  22)2  [z->oo  (7r‘  —  22)“ 


C0S(2)  =  0, 


C0S(2)  =  2C0S(2), 


where  in  the  last  two  lines  we  formally  separate  the  oscillating 
(and  bounded)  cosine  factors  in  order  to  extract  th e  functional 
form  of  the  surviving  terms.  Consequently, 

Fz(uj)  «  2  (1  +  cos(W)) ,  wr>  1,  (D18) 


and  the  integral  is  well  approximated  by 


2  (1  +  cos(wr))  duj 

2+  sin(a;cT)  \ 
wcr  ) 


«  2 U!c,  (VCT  1. 


(D19) 

(D20) 

(D21) 
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